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ABSTRACT 


During July 2007, fields from both the NCEP GFS final analyses and 
ECMWF model of OW and SF analyses suggests the development of two distinct 
areas with sub-tropical intrusion from remnants of a decaying baroclinic system 
in the WNP. 

This analysis of the formation of Usagi points to sub-tropical intrusion of a 
strong lower-tropospheric baroclinic system undergoing decay as potential 
seedlings for typhoon formation in areas of high sea surface temperatures, weak 
low-level vertical wind shear, and persistent convection. As the PV anomalies is 
stretched and detached from the baroclinic source region, it is wrapped around a 
strong tropospheric anticyclone in the subtropics. This constitutes a different 
type of baroclinic initiation process than has been previously identified in Atlantic 
cyclone formation events associated with TT, which are induced by upper-level 
troughs. 

The area of high values of OW at the tip of a PV streamer favors 
sustained deep convection, which will enhance the low-level vorticity and moisten 
the mid-level thereby producing high values of SF. The area of strong vorticity at 
the tip of the second PV streamer possessed both high OW and SF, favoring 
deep convection and cyclonic vortex tube stretching that appeared to culminate 
in an enhancement of lower tropospheric cyclonic vorticity. 

Although this analysis was originally motivated by initial analyses 
suggesting that ex-hurricane Cosme underwent a direct vorticity interaction with 
the second PV streamer, our revised hypothesis on the role of Cosme is that it 
may have enlarged the wave pouch and helped preserve the mid-tropospheric 
circulation from hostile outside influences. In this sense, the Cosme wave may 
have played an important indirect role in the formation of Usagi. This hypothesis 
requires further investigation. 
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I. INTRODUCTION 


Improved forecast of tropical cyclogenesis is vital to US military interest in 
the western North Pacific (WNP) and other maritime activities in this region. A 
weather phenomenon such as a tropical cyclone (TC) can have dangerous 
winds, flood-causing heavy rain, and destructive storm surge over large areas so 
the degree of impact can be high. An approaching Typhoon can curtail the 
civilian economy, military operations, and all routine day-to-day living along 
several hundred miles of coastline and inland areas. The duration of the impact 
can be several days to several months depending on the damage. 

During 22 -28 July 2007 two distinct areas of disturbed weather 
associated with a subtropical intrusion (SI) emerged out of remnants of a 
decaying baroclinic system in the WNP. In this thesis, key events that led to the 
non-development of the initial invest area and the development of Typhoon Usagi 
(5W) are investigated. This work will explore the usefulness of satellite data and 
model analysis to include the first testing in the WNP of Okubo-Weiss (OW) 
criterion and saturation fraction (SF), which are parameters that diagnose factors 
hypothesized to have aided in the formation of Typhoon Usagi. 

The genesis of Usagi was initially analyzed using the GFS model and the 
analysis of relative vorticity at 850 hPa with the VORTRACK vorticity tracking 
system. The initial hypothesis was that Usagi was a tropical transition (TT) of a 
strong equatorial cyclone (SEC) that moved southwest as a decaying mid¬ 
latitude baroclinic system. It was then hypothesized that Usagi intensified due to 
vortex merger and axisymmetrization between the tail-end of the decaying 
baroclinic system and remnants of Hurricane Cosme (5E) that formed in the 
eastern North Pacific and moved westward across the dateline. A more detailed 
discussion of this possible sequence will follow in Chapter III. Based on high- 
resolution model analysis, this initial hypothesis was revised to include the 
circulation of Cosme as a factor that widened a “pouch” of vorticity anomalies of 


1 



Usagi. Some preliminary tests of this hypothesis are presented using the model 
analyses. The high time-resolution analysis will demonstrate no direct vorticity 
connection between the two systems; and satellite imagery reveals a persistent 
area of clearing between the two systems that is consistent with a region of 
subsidence. This evidence suggests that the wave-like remnants of Cosme are 
too far away from the baroclinic system (its closest proximity was approximately 
1000 km on 0600 UTC 26 July) to be a deciding influence in the formation of 
Usagi. This study suggests an alternative hypothesis that the enlarged wave 
pouch created a favorable condition for the genesis of Usagi, especially at 600 
hPa. 

The tip of the potential vorticity (PV) streamer, which is dominated by 
curvature vorticity, is believed to be the key factor leading to the formation of 
Typhoon Usagi. A natural question is whether or not Usagi would have 
eventually developed without Cosme in widening the vorticity pouch. These 
questions are beyond the scope of this paper and should be the focus of further 
studies. 

This study will highlight the difference between Atlantic and WNP TT 
cases and suggest why the TT model as defined in the literature for the Atlantic 
sector fails to give useful forecast rules for the WNP. This case offers a new 
example of a genesis sequence from an extratropical system in the WNP. While 
the frequency of this weather sequence is believed to be relatively low, the 
formation of any typhoon can have severe impact. 

The series on TC genesis described in Dunkerton, Montgomery and Wang 
(2008; hereafter DMW08) examined the TC genesis pathway in association with 
easterly waves and instability of the Inter-tropical Convergence Zone in the 
Atlantic and eastern Pacific sectors. Their three-year analysis of developing 
cases in these sectors revealed a common occurrence of the creation of a closed 
gyre or “pouch” that is protected to some degree from hostile environmental 
influences. 
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This analysis of the formation of Usagi points to SI of a decaying baroclinic 
systems as potential seedlings for typhoon formation in areas of high sea-surface 
temperatures, weak low-level vertical wind shear, and persistent convection. In 
the case of Usagi, a large-scale anticyclone alone in the west side of the 
baroclinic zone is apparent. This work will use MTSAT satellite, PV, OW, SF, 
and TRMM retrieval analyses to show that the stretched PV anomalies that get 
detached from their baroclinic source region and subsequently wrapped around 
the large anticyclone is a different type of baroclinic initiation mechanism than the 
TT Atlantic cyclone formation model, which are induced by upper-level troughs. 
Thus, the cyclogenesis of Usagi in the WNP is not consistent with TT cases in 
the Gulf of Mexico or Atlantic region. 

Prior to the formation of Usagi, an initial disturbance was labeled as an 
invest region for possible tropical cyclogenesis formation. However, this non¬ 
developing system exhibited less organization and weaker vorticity than the 
developing (Usagi) case. The failure of the initial SI to develop may occur due to 
the cyclone translating over cool waters, an area of strong vertical wind shear, a 
lack of moisture, or insufficient cyclonic vorticity, as the cyclone becomes 
detached from the main vorticity reservoir in the decaying baroclinic system. 

A. ATLANTIC TT AND WNP CYCLONE FORMATION 

According to Davis and Bosart (2004; hereafter DB04), TT refers to the 
formation of a warm core tropical cyclone from a cold-core disturbance, usually of 
extratropical origin. Some of the largest forecast errors of track and intensity 
occur with these events. Their proximity to the North American coast produces 
significant adverse consequences for poor forecasts, namely, minimal time for 
coastal warnings and inland flood preparations. DB04 define two TT scenarios 
based on the amplitude and structure of the precursor seedling disturbance: SEC 
and weak equatorial cyclone (WEC). In SEC scenarios, the extratropical 
cyclogenesis produces a surface cyclone capable of wind-induced surface heat 
exchange (WISHE). In contrast, WEC situations have the baroclinic cyclone as 
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the organizing agent for convection. This convection must then undergo self¬ 
organization to create a disturbance capable of self-amplification. DB04 suggest 
no clear threshold exist between SEC and WEC cases, since these cases 
represent points on a spectrum. Once an appropriately strong surface vortex is 
formed, no clear distinction was found in the ensuing tropical cyclone 
intensification in either SEC or WEC situations. 

In these stronger extra-tropical cases, baroclinic and diabatic processes 
combine to produce an intense mesoscale cyclone, which then transitions to a 
tropical cyclone (Davis and Bosart 2003, hereafter DB03) as the WISHE process 
develops a low-level, warm-core circulation center. DB04 showed that SEC 
cases have a more consistent and repeatable evolution than WEC cases, and 
strong tropical cyclones tend to form most often with SEC scenarios. 

In contrast, the baroclinic WEC situations act as the organizing agent for 
convection. This convection must then undergo self-organization to create a 
disturbance capable of self-amplification. DB04 suggest no clear threshold exist 
between SEC and WEC cases, since these cases represent points on a 
spectrum. Once an appropriately strong surface vortex is formed, no clear 
distinction was found in the ensuing tropical cyclone intensification in either SEC 
or WEC situations. 

The TT cases discussed above reside within a broader continuum of 
marine cyclogenesis that comprise a range from baroclinic cyclones to 
hurricanes initiated from weak extratropical systems (DB04). Although storm- 
induced surface heat and moisture fluxes play a key role in TT, the key role in TT 
seems to be determined by baroclinicity as long as the circulation is over a 
sufficiently warm ocean. However, forecasting the TT of various baroclinic 
systems in different areas such that of this WNP case is rather challenging. 

In SEC instances, the low-latitude frontal cyclone has intensity sufficient to 
begin the WISHE process. DB03 showed diabatic processes serve to eliminate 
the vertical wind shear present over the cold-core low. The decrease in shear is 
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caused by both the upper-tropospheric outflow from convection, and from the 
diabatic redistribution of PV. Both of these processes help to homogenize the 
horizontal gradients of PV directly above the low-level circulation. As a result, an 
equilibrated cyclone resembling an occluded system in a sub-synoptic “cocoon” 
of weak shear is generated. This Usagi case takes place within this favorable 
environment, and eventually a tropical cyclone forms. 

The occlusion before TT is different from the classical model of occlusion 
in which surface cyclone moves poleward toward colder air beneath an upper- 
tropospheric jet. In contrast, the jet in the TT cases is rearranged via diabatic 
processes and results in an apparent migration of the surface cyclone toward 
warmer tropospheric air (DB04). 

The rainfall and cloud signatures just before TT began in the four SEC 
cases discussed in DB04 are reproduced here in Figure 1. In each situation, a 
marked asymmetry exists in the rainfall, with a proclivity for a “bent back” frontal 
structure, and heavy rainfall on the west side of the still extratropical surface low. 
This structure was observed prior to some of the stronger hurricanes generated 
from the TT process. The pre-hurricane Florence structure in Figure la most 
resembles the pre-typhoon Usagi seedling. 
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85 GHz Polarized Corrected Temperature (PCT) (K) 


Figure 1. (a) Visible satellite and Tropical Rainfall Measuring Mission (TRMM) 
Microwave Imager (TMI) 85-GHz polarized corrected temperature (PCT) imagery 
at 1145 and 1209 UTC 10 Sep 2000, respectively; (b) (upper right), as in (a) but 
for visible and TMI 85-GHz PCT at 1245 and 1335 UTC 15 Oct 2000; (c) IR and 
Special Sensor Microwave/Imager 85-GHz PCT at 0015 and 0219 UTC 10 Sep 
2002; and (d) IR and TMI 85-GHz PCT at 2045 and 2132 UTC 11 Oct 2001. In all 
panels, the “L” indicates the position of surface low. Images obtained courtesy of 
Naval Research Laboratory: 14 February 2008. 

In WEC cases, precursor disturbances can be weak baroclinic waves, or 
mid-tropospheric vortices that simply have insufficient amplitude to create a 
surface cyclone capable of amplifying by WISHE without invoking an 
intermediate process to enhance mesoscale vorticity. The mechanism 
responsible for the warm-core transition in WEC cases is unclear. The 
production of cyclonic vorticity centers on scales of 10-20 km, so-called rotating 
hot towers (Hendricks et al. 2004, Montgomery et al. 2006), may be the preferred 
means by which the convective scale influences the development of the primary, 
larger-scale vortex. Mid-tropospheric convectively generated vortices have also 
been observed to initiate tropical cyclogenesis in lower latitudes (e.g., Simpson et 
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al. 1997). DB04 argue that it is through understanding the initiation of tropical 
cyclogenesis by mesoscale vortices at higher latitudes that the link with tropical 
cyclone formation in the deep tropics can be made. Just as there is a continuum 
between the weak and strong baroclinic precursors, there is a range of 
realizations between mesoscale vortices and weak baroclinic systems. 
Numerical simulations of Diana (1984) (Powers and Davis 2002) showed that the 
path to tropical cyclogenesis required a lower-middle-tropospheric vortex that 
formed within convection initiated by an extratropical precursor. The 
distinguishing character between cases, the mid-tropospheric vortex existed for 
more than two days before tropical cyclogenesis began. In such cases, the 
vortex must organize new convection through its interaction with vertical shear. 
When a weak extratropical cyclone is maintained, the organization of convection 
occurs through a superposition of vortex and cyclone-induced ascent. In the 
cases studied, there were multiple lower tropospheric vortices formed within a 
mesoscale ascent region, and the coalescence and growth of these vortices 
formed the nascent tropical storm (TS) (Hendricks et al. 2004). Considerable 
research remains to understand how convection organizes in systems with weak 
precursors. Therefore, it is difficult to derive a set of forecast rules. However, it is 
apparent that PV “debris” extruded from the mid-latitude jet is common over the 
warm oceans of the subtropical Atlantic, even as far south as 15°N on occasion. 
Overall then, favorable conditions for WEC cases of TT involve mid-upper- 
tropospheric cyclone PV anomalies encountering lower-tropospheric baroclinity 
(and, hence, vertical shear). In the cases that develop into tropical cyclones, the 
shear remains either roughly 10 m s” 1 or less, or it is reduced to such values 
following convection organization (Powers and Davis 2002). DB04 noted that at 
lower latitudes, convection has also been observed to organize when upper- 
tropospheric disturbances approach easterly waves, wherein systematic 
baroclinity is hard to identify. There are almost certainly other factors involved in 
WEC cases of TT, but research has yet to fully clarify them. Although in the 
Gabrielle (2001), the period is relatively short (about 6 hours); the depression 
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undergoes a full warm-core transformation (not shown). It is believed that the 
ongoing convection is somehow responsible for this decrease in shear. 

One of the critical forecast challenges, then, is to anticipate a favorable 
environmental modification, since it usually occurs on the timescales of less than 
one day. DB04 proposed a primary empirical forecast rule to aid in the forecast 
challenge of predicting SEC TT cases: “The precursor cyclone must occlude and 
remain over warm water (>26°C) for at least a day after the occlusion.” DB03 
hypothesized that the initially strongly sheared environment, which is hostile to 
tropical cyclone formation, is significantly modified as the incipient cyclone 
progresses. After about one day, the vertical shear is markedly decreased, 
which facilitates the creation of a finite amplitude vortex in the lower troposphere 
that is able to enhance the moist entropy fluxes from the ocean surface. 

Hurricane Diana provided a useful illustration of this process. Diana was a 
modest baroclinic cyclogenesis event that stimulated organized convection, 
which led to the coherent vortices in the lower troposphere that formed the seed 
for tropical cyclogenesis (DB03, Hendricks et al. 2004). Another prime example 
of TT was in the case of Hurricane Michael (2001) (not shown). The pre-Michael 
cyclone was a robust baroclinic cyclone that was able to generate a finite 
amplitude mesoscale surface cyclone so the transition to a warm-core system 
could begin almost as soon as the vertical shear weakened. 

Failure of TT may occur because of two reasons. Transition will fail if the 
occluding cyclone translates over cool water before TT can occur (thus, the “1- 
day rule” above). Secondly, transition will fail if the primary cyclone is prevented 
from occluding due to an increase in vertical wind shear over its center from short 
waves approaching the surface cyclone (DB04). 

B. THE MARSUPIAL PARADIGM 

Summarizing the DMW08 theory, the relevant theoretical paradigm is the 
formation of a closed proto-vortex or “embryo” (Emanuel, 2007; Montgomery, 
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2006) within the wave in the lower troposphere. Subsequent eddies that are 
shed from the vortex, become energized further by deep convective heating to 
become detached from the baroclinic wave. In the following section, we will refer 
to this as the marsupial paradigm. The marsupial paradigm is a rudimentary 
understanding of the Lagrangian flow and its dynamical and thermodynamical 
properties in the neighborhood of the developing storm. 

DMW08 proposed three main hypotheses for TC formation that 
synthesized existing knowledge and advanced new insight on the flow kinematics 
and dynamics, moist thermodynamics, and wave/vortex interactions during the 
transformation into a tropical depression-strength vortex. First, the proto-vortex 
cyclonic eddies instrumental in TC formation are hypothesized to be intimately 
associated with the parent wave’s critical latitude in the lower-troposphere. The 
critical level/layer is generally a surface in three dimensions, but the word level 
has been used in this study to represent one of these dimensions, e.g., latitude. 
This structure forms a critical layer because of the wave’s finite amplitude 
interaction with its own critical latitude. The critical layer is a region of cyclonic 
rotation and weak strain/shearing deformation in which synoptic and mesoscale 
anomalies move together and amplify on a nearly zero relative mean flow. This 
multi scale interaction provides a dynamical pathway to bottom-up development 
of the proto vortex. The second hypothesis is that the critical layer of the parent 
wave provides a set of closed material contours inside of which air is repeatedly 
moistened by convection. The wave is protected to some degree from lateral 
intrusion of dry air impinging vertical shear, and able to keep in pace with the 
parent wave until the proto-vortex has strengthened into a self-maintaining entity. 
The third hypothesis is that the parent wave is maintained and possibly 
enhanced by diabatically amplified eddies within the wave. This process is 
favored in regions of small intrinsic phase speed. This hypothesis, agrees with 
the common observation that the tropical wave is considerably weakened or 
eliminated when the diabatic vortex leaves the pouch, having acquired its own 
identity and propagation characteristics independent of the parent wave. 
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The ex-hurricane Cosme disturbance is an illustrative example in which a 
westward propagating wave carried its proto-vortex over time. Once the 
depression ensued, the strengthened and enlarged vortex (Usagi) acquired a 
poleward component of propagation while its zonal component begins to slow 
relative to the phase speed of the parent wave. 

C. ELEMENTS OF THE GENESIS SEQUENCE 
1. Dynamical Fields - OW Parameter 

The relative vorticity ) and OW fields at 600 hPa and 850 hPa will be 
examined to single out curvature vorticity compared to shear vorticity. The large 
positive values of OW parameter suggest regions where curvature vorticity is 
dominant. Perhaps the most important conclusion to draw from the OW 
parameter inside the translating gyre is that the flow within it is relatively free of 
strain and shear (McWilliams 1984), favoring nominally the formation of smaller 
proto-vortical structures within. 

The moving frame of reference is optimal to display the Lagrangian flow 
and the evolution of the storm vortex. This enables one to estimate the 
Lagrangian parcel motions that are affecting the dynamics and thermodynamics 
of the transformation process. 

The OW parameter is defined as 

ow = C~s 2 ( 1 . 1 ) 

where 

du dv\ (dv du 

-+- 1 - 

dx dy) dy 

is the square of the horizontal strain deformation and 



( 1 . 2 ) 
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is the relative vertical vorticity. 

Negative values of the Okubo-Weiss parameter for the simulation period 
are regions that are strain dominated, whereas positive values indicate strain- 
free rotationally dominant regions. The vortex cores are characterized by 
strongly positive values of OW (Weiss, Pasquero, and Provenzale 2002). In the 
following, we identify the vortex cores as those circular regions with positive OW. 

2. Thermodynamical Fields - SF 

The saturation fraction of the lower troposphere, defined as the ratio of 
precipitable water to saturated precipitable water, from case studies concluded 
by David Raymond (2002), the saturation fraction appears to be the primary 
control on precipitation in the ECMWF model and the surface moisture flux is 
acting in a secondary fashion via its effect on the SF. In this thesis, we use 600 
hPa SF as a benchmark to focus on the lower tropospheric moisture anomalies. 

D. OBJECTIVES OF CURRENT STUDY 

1. To present the facts that led to the development of Usagi; examine 
model analysis data, and satellite imagery in observing the conditions leading to 
WNP TC formation 

2. To explore the usefulness high-resolution model analyses of OW and 
SF in order to assess these hypothesized factors in the formation of Typhoon 
Usagi (5W). 

3. To propose an explanation on how the remnants of Hurricane Cosme 
(6E) created a favorable condition for the genesis of Usagi. 

More generally, we wish to increase the forecast capability of the Joint 
Typhoon Warning Center (JTWC) and we believe that further investigation and 
qualitative documentation of the various pathways to tropical cyclone formation 



will prove useful to the forecasting community as well. This includes identifying 
synoptic patterns favorable for tropical cyclone formation and providing guidance 
on mesoscale conditions that may enhance or inhibit tropical cyclone formation in 
the WNP (Non-developer and developer cases). 

The organization of the thesis is as follows. In Chapter II, an overview of 
the satellite and model data is presented, and this is followed by discussion of 
the data processing used in the thesis. In Chapter III, the case studies used in 
the thesis are presented. Finally, the summary of results, conclusions, and 
recommendations for future research are contained in Chapter IV. 

E. PRE-STORM SITUATION 

1. Hurricane Cosme (6E) 

On 14 July 2007, the National Hurricane Center (NHC) identified Tropical 
Depression (TD) (6E) as an area of disturbed weather about halfway between 
Mexico and the Hawaiian Islands. The area of disturbed weather acquired a 
surface circulation and sufficiently deep convection for the NHC to designate it as 
a TD, what then became TS Cosme. 

On 24 July, the JTWC was still tracking the weak remnants of Cosme 810 
nautical miles east of Wake Island. Animated infrared satellite imagery depicted 
improved consolidation with persistent deep convection near the partially 
exposed low-level circulation center. 

2. Baroclinic Extratropical Low 

This extratropical low moving southwest out of the mid-latitudes into the 
tropics had a low-level circulation that was already near the surface and rather 
strong. Thus, it did not take much to transition this cold-core low into a tropical 
system given that it was moving southward into an area with warm sea-surface 
temperatures (SSTs), high ocean heat content, low vertical wind shear, and the 
approaching remnants of Cosme. 
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On 22 July, an area of deep convection emerged out of this system but 
never developed. However, on 28 July, JTWC issued a Tropical Cyclone 
Formation Alert (TCFA) from the same extratropical low, and then designated the 
system a TD later that day. 

3. Tropical System Non-developer (Gyre) 

As mentioned earlier, two noted systems of deep convection developed 
from remnants of a decaying baroclinic extratropical low. The initial system on 
22 July gets detached but never develops and the second identified system on 
28 July eventually became Typhoon Usagi. To look at the two systems, a simple 
strategy is to compare the precursors of the two tropical system formations 

4. Tropical Cyclogenesis of Typhoon Usagi (5W) 

According to JTWC (2007), on 1200 UTC 28 July the system quickly 
strengthened as a TD as it approached 260 nautical miles northeast of Saipan 
and had tracked westward at 14 knots over the past 6 hours with maximum 
sustained winds at 25 knots gusting to 35 knots. At 0000 UTC 29 July, JTWC 
upgraded it to TS strength. 

One of the most interesting aspects of tropical cyclogenesis that is not 
understood, especially in the western North Pacific (WNP) sector, is SI or TT 
events. According to DB04, the development of almost half (26 out of 57) of all 
Atlantic tropical cyclones from 2000 to 2003 relied on an extratropical 
disturbance. Majorities of these systems were of a baroclinic origin and were 
initially cold-core systems. A kinematic and thermodynamic (cold core to warm 
core) transformation of these disturbances was required. DB04 consider this 
process to be the definition of TT, and contrast it to extratropical transition (ET), 
in which a TC transitions into an extratropical cyclone. However, this has not 
been considered as a dominant formation type in the WNP. The thesis will 
examine the formation of Usagi relationship to the TT conceptual model. 
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II. DATA AND METHODOLOGY 


The General Meteorological Package (GEMPAK) Analysis and Rendering 
Program (GARP) was used to generate plots of fields from the National Center 
for Environmental Prediction (NCEP) Global Forecast System (GFS) 30-day 
global analysis. Additionally, the Navy Operational Global Atmospheric 
Prediction System (NOGAPS) model analysis and the United Kingdom 
Meteorological Office (UKMET) model were also examined. All of the model data 
was examined from 0000 UTC 22 July through 1800 UTC 8 August 2007. 

At a later stage of this research project, the NCEP GFS global analysis 
Global Data Assimilation System (GDAS) final analysis (FNL) was added in order 
to examine PV and streamlines. To calculate the OW and SF parameters, the 
European Centre for Medium-Range Weather Forecasts (ECMWF) model data 
was used. 

Vis5D supplements the GARP plots for more detailed analysis. Vis5D is a 
software system that can be used to visualize both gridded data and irregularly 
located data. Sources for these data can come from numerical weather models, 
surface observations and other similar sources. Vis5D can work on data in the 
form of a five-dimensional rectangle. That is, the data are real numbers at each 
point of a grid, which spans three space dimensions, one time dimension and a 
dimension for enumerating multiple physical variables. Of course, Vis5D works 
perfectly well on data sets with only one variable, one time step (i.e. no time 
dynamics) or one vertical level. However, the data grids should have at least two 
rows and columns. The Vis5D can also work with irregularly spaced data that 
are stored as records. Each record contains a geographic location, a time, and a 
set of variables that can contain either character or numerical data. A major 
advantage of Vis5D is that it readily supports intercomparison of multiple data 
sets. Data can be overlaid in the same 3-D display and/or viewed side-by-side 
spread sheet style. Data sets may be overlaid and aligned in space and time. In 
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the spreadsheet style, multiple displays can be linked. Once linked, the time 
steps from all data sets are merged and the controls of the linked displays are 
synchronized. 

The University of Wisconsin - Cooperative Institute for Meteorological 
Satellite Studies (UW-CIMSS) and the National Oceanic and Atmospheric 
Administration - Hurricane Research Division (NOAA- HRD) products were also 
used to supplement the satellite data used here. This imagery was used as an 
overlay on the previous GARP generated GFS 0000-hour analysis fields of 850 
hPa winds (knots) and 200 hPa winds (knots), and 850 hPa absolute vorticity 

(with values greater than 10 5s ’shaded) from 0000 UTC 22 July through 1800 
UTC 28 July. The imagery overlay helps contrast where the convective activity 
is located versus the analyzed generated absolute vorticity. 

A smaller area within in the synoptic region was chosen and ingested 
inside the Vis5D software in order to simulate trajectories where the convection 
and vorticity was initiated. The streamlines were superimposed in the wave¬ 
centric perspective. The Vis5D analysis region was also used to obtain a sense 
of the area of deep convection in comparing the developing and non-developing 
circulation. Plots of the low-level convergence vs. cold cloud-top temperatures 
were constructed (Vis5D, 2005). 

The Air Force Weather Agency’s (AFWA’s) Interactive Grid Analysis and 
Display System (IGrADS) visualization tool was also used to display the detailed 
graphics. The Center for Ocean-Land-Atmosphere Studies (COLA) IGrADS is 
available over the AFWA’s webpage, Joint Air Force and Army Weather 
Information Network (JAAWIN) (https://weather.afwa.af.mil) and has been 
operational since early 2002. The interface enables forecasters to create 
varieties of visualized or alphanumeric products designed to meet their 
customers’ needs. These products range from meteograms, skew-Ts, vertical 
cross-sections, forecast maps, and seven different types of alphanumeric 
meteorological output. Forecasters have the option of choosing from among 
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eight different numerical forecast models. The most extensive product lines are 
available for AFWA’s MM5. IGrADS users can choose any display region they 
desire, can animate forecast data, and can bookmark frequently needed charts, 
and then recall them with current model data whenever required (COLA, 2004). 

A. SATELLITE DATA PROCESSING 

1. MTSAT Infrared-Red Imagery 

Multi-Functional Transport Satellite (MTSAT) series fulfill two functions: a 
meteorological function by the Japan Meteorological Agency and an aviation 
control function by the Civil Aviation Bureau of the Ministry of Land, Infrastructure 
and Transport. The MTSAT series succeed the Geostationary Meteorological 
Satellite (GMS) series as the next generation satellite series covering the East 
Asia and the Western Pacific regions. Satellite images are available over the 
WNP region and over each of the tropical western Pacific sites. The Cloud-top 
temperatures defined from MTSAT imagery will be used with thresholds to 
distinguish between warmer lower level clouds (from -20°C to less than -70°C) 
and colder upper-level cloud tops. 

The raw data are down-inked to the SeaSpace Corporation receiving 
station at the University of Hawaii. GMS-5 was launched in March 1995 into a 
geostationary orbit above 0° N, 140°E. It carries a Visible and Infrared Spin Scan 
Radiometer (VISSR) with four channels, one in the visible part of the spectrum 
and three in the infrared. The ground resolution at the sub satellite point is 1.25 
km in the visible channel and 5 km in the infrared. A set of matching overlay files 
showing a latitude and longitude grid and the positions of the coastline are 
available for the WNP 5km resolution spatial coverage. The overlay files consist 
of an image of value one everywhere, except at the positions of the features 
(coastlines) in which case the pixel value is zero. Thus, the overlays can be 
written into the image by multiplying them together; the image data remain 
unchanged and the overlain pixels become zero. 
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2. Tropical Rainfall Measuring Mission (TRMM) 


Accurate measurement of the spatial and temporal variation of tropical 
rainfall around the globe remains one of the critical unsolved problems of tropical 
meteorology. TRMM, with its broad sampling footprint between 35°N and 35°S, 
is providing some of the first detailed and comprehensive dataset on the four 
dimensional distribution of rainfall and latent heating over vastly under sampled 
oceanic and tropical continental regions. Combined with concurrent 
measurement of the atmosphere's radiation budget, estimates of the total 
diabatic heating are being realized for the first time ever on a global scale. 
TRMM will fill many gaps in our understanding of rainfall properties and their 
variation. These includes frequency distributions of rainfall intensity and area 
coverage; the partitioning of rainfall into convective and Stratiform categories; the 
vertical distribution of hydrometeors (including the structure and intensity of the 
stratiform region bright band), and variation of the timing of heaviest rainfall. 
These are particularly relevant in the documentation of nocturnal intensification of 
large mesoscale convective systems over the oceans, and diurnal intensification 
of orographically and sea-breezed forced systems over land. TRMM will enable 
mapping of larger time and space variations of rainfall in quasi-periodic 
circulation anomalies, such as the Madden-Julian Oscillation (MJO) in the WNP 
and El Nino Southern Oscillation (ENSO) over the broader Pacific basin. 
Furthermore, the critical onset of large annual circulation regimes, such as the 
Asian summer monsoon, can be more thoroughly studied (TRMM 2008) 

B. NUMERICAL MODEL FIELDS 

This case study utilizes a variety of numerical model fields from 
operational global models, and derived products to analyze the evolution of the 
non-developer and the Sis that became Typhoon Usagi during 22-28 July 2007. 
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1 . 


GFS 


The NCEP GFS global analysis GDAS FNL, used as initial conditions for 
the GFS, will be used primarily. The GFS 0000, 0600, 1200, and 1800 UTC 
analysis fields were used to provide a 6-hour time snapshots. Only the 00-hour 
analyses will be used to avoid model forecast bias from entering the forecast 
plots. The 850 hPa and 200 hPa streamlines will be used to illustrate both the 
synoptic flow patterns at their respective levels. The 850 hPa absolute vorticity 
(with values greater than 10'V 1 shaded) will be used to correlate with active 
areas of deep tropical convection. Finally, the distribution of PV within the levels 
of 850-500 hPa will be used to illustrate low-level circulation. The potential 
vorticity will be displayed in potential vorticity units (PVU), defined as 10' 6 m 2 s _1 K 
kg' 1 following Hoskins et al. (1985). The thresholds used in GARP are at a 
minimum of 0.3 and a max of 0.7 x 10' 6 m 2 s _1 K kg" 1 with intervals of 0.1 X 10" 6 
m 2 s _1 K kg' 1 . 

2. NOGAPS and UKMET 

The NOGAPS and UKMET model analysis fields provided by Naval 
Research Laboratory-Monterey (NRL-MRY) are used in this study at 0000 UTC 
and 1200 UTC to examine the baroclinic system, TC structure and dynamics as 
compared to IR, WV, and when available, microwave imagery. Because the 
0600 UTC and 1800 UTC analyses were not available, the 6-h forecasts from the 
0000 UTC and 1200 UTC model integrations respectively were substituted to 
provide a 6- hour time resolution. The 850 hPa and 200 hPa streamlines will be 
used to illustrate both the synoptic flow patterns at their respective levels. The 
850 hPa absolute vorticity (with values greater than 10'V 1 shaded) will be used 
to correlate with active areas of deep tropical convection through the whole TT 
period. 
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3. 


ECMWF 


The ECMWF forecasting system consists of five components: a general 
circulation model, an ocean wave model, and data assimilation. In 1998, a 
seasonal forecasting system started to operate and in 2002, a monthly 
forecasting system was introduced (ECMWF, 2007). 

C. LIMITATIONS 

The limitations of the data sets used are model resolution, data gaps, and 
the time of satellite imagery. The GFS model horizontal resolution is spectral 
triangular 254 (T254) with a Gaussian grid of 768 x 384, which is roughly 
equivalent to 0.5° x 0.5° latitude/longitude. The vertical domain is divided into 64 
layers with enhanced resolution near the bottom and the top. For a system with 
a surface pressure of 1000 hPa, the lowest atmospheric level in the model would 
be at a pressure of about 997.3 hPa, 15 levels would be below 800 hPa, 24 
levels would be above 100 hPa and the top level would be at about 0.27 hPa 
(NCEP, 2004). By contrast, the NOGAPS model is 0.75 degree horizontal 
resolution (spectral resolution = T159) and 24 levels in the vertical, from the 
Earth surface to 50 km. The model runs operationally on 12 processors with a 
sustained performance of five Gflops (5><10 9 FLoating point Operations Per 
Second) (NRL, 2003). 
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III. ANALYSIS AND DISCUSSION 


The scope of this study will be from the inception of Cosme until the 
weakening of Usagi with focus on 22-28 July 2007. The spatial domain is from 
longitude 130°E to 179°W and latitude 30°N to 10°N. The study will demonstrate 
that the tip of the PV streamer, which is dominated by curvature vorticity, is the 
key factor leading to the TC formation. 

The main objective of this study is to show the comparison between the 
non-developing and developing (Usagi) systems. A simple strategy is to 
compare the precursors of the two tropical systems by looking at the large-scale 
tropospheric environment of Usagi, lower tropospheric interactions within the 
non-developer (gyre) and developer (Usagi), by employing the OW and SF 
parameters. 

The GFS initial model analysis was examined, using the VORTRACK 
system analysis (Figure 2-4) that identifies vorticity centers and tracks them. 
Examination of this model data led us to the initial hypothesis that the Usagi 
system resulted from the vortex merger of two distinct convective systems 
(Cosme and Usagi) of different origin. The interaction between the two systems 
was originally believed to have led to a warm-core system that eventually 
became typhoon (Usagi). The VORTRACK analyses suggest that as the two 
vortices approached each other, they start to undergo axisymmetrization and 
filamentation (Figure. 2-4). Initial analysis using the VORTRACK model led us 
believe that Usagi formed by the ingestion of like-sign PV anomalies 
(Montgomery and Enagonio 1998). 


21 




Non-developer 
= 5 = 


Usagi 


Cosfjie wave 


170 W 


15N 

ION 


40N 

35N 

30 N 

25N 


2007072612+000 


Baroclinic System 


Figure 2. VORTRACK analysis of Usagi and Cosme interaction on 1200 UTC 26 
July 2007 48 hours before Usagi was named as TD 05W. Shadings are relative 
vorticity (unit: 1X1 O' 5 ), winds @ 850 hPa plotted if relative vorticity is > 0, and 
black dots are track (12 hour steps) where relative vorticity is 1.5 (unit 1.5 X 10' 5 ). 
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Figure 3. VORTRACK analysis of Usagi and Cosme interaction on 0000 UTC 27 
July 2007 36 hours before Usagi was named as TD 05W. Shadings are relative 
vorticity (unit: 1X1 O' 5 ), winds @ 850 hPa plotted if relative vorticity is > 0, and 
black dots are track (12 hour steps) where relative vorticity is 1.5 (unit 1.5 X 10' 5 ). 
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Figure 4. VORTRACK analysis of Usagi and Cosme interaction on 1200 UTC 27 
July 2007 24 hours before Usagi was named as TD 05W. Shadings are relative 
vorticity (unit: 1X1 O' 5 ), winds @ 850 hPa plotted if relative vorticity is > 0, and 
black dots are track (12 hour steps) where relative vorticity is 1.5 (unit 1.5 X 10' 5 ). 

When the GFS FNL and ECMWF model analysis became available, a 
closer investigation of the high resolution analysis data using OW and SF 
parameters lead to a revision to the initial hypothesis and an explanation of the 
events leading to the Usagi TC formation as being fundamentally different from 
vortex merger, axisymmetrization, and filamentation. This study will document 
the difference of TT in the WNP. 

A. OVERVIEW OF INITIATING BAROCLINIC SYSTEM 

1. The Large-Scale Lower Tropospheric Environment of Usagi 

From 22-28 July 2007, two distinct areas of disturbed weather with SI 
emerged out of remnants of a decaying baroclinic system in the WNP. The initial 
system on 22 July becomes detached but never develops and the second 
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identified system on 25 July eventually became Typhoon Usagi. In the non¬ 
developing system, an area of weaker vorticity and a Tropical Upper 
Tropospheric Trough (TUTT) damped the convection. In the Usagi case, a 
strong line of low-level confluence and convection occurs along the line of 
confluent region. This region of convection will continue to intensify until the 
whole system develops as an organized center of deep convection. 

The PV streamers are detached from the baroclinic system as the upper- 
level low is occluded. This formation precursor is still cold core since it is a 
reflection of an upper-level trough axis. As time progresses, this PV anomaly 
gets detached and is wrapped into the anticyclone. This baroclinic type of 
initiation process may not be seen in the Atlantic cyclone signatures from past 
research studies (DB04), which are induced by upper-level troughs. In this case, 
a large-scale anticyclone is to the west of the baroclinic system. In the lower 
levels, the vorticity anomaly becomes organized with enhanced low-level 
convergence from the existing instability from the low-level. This would be 
consistent with the bottom-up formation of a frontal system stagnating in the 
subtropics and transforming to warm-core characteristics. 

The synoptic features indicate areas of strong enhanced convection had 
been local and mainly associated with tropical cyclones north of the Philippines. 
From 16 July to 1 August, the MJO signal remained weak (Figure 5). During 
these weeks, convective anomalies had generally been weak across the deep 
tropics, with slightly above average conditions across the eastern Indian Ocean 
and western Maritime Continent. The outgoing long-wave radiation anomalies 
during late July and early August confirmed the wet conditions over the Maritime 
Continent. 

From mid-June into early August, weak to moderate MJO activity was 

observed as velocity potential anomalies propagated eastwards (Figure 6). 

Since mid-June, weak-moderate MJO activity had been observed as regions of 

suppressed and enhanced convection had shifted eastward from the Indian 

Ocean into the far western Pacific. Convection had become near average across 
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much of the equatorial tropics with dry conditions only evident just west of the 
date line. The 850 hPa zonal wind easterly anomalies increased during mid-July 
over the Maritime continent and western Pacific and shifted eastward during mid- 
late July. The MJO index is defined by Wheeler and Hendon (2004) based on 
the first two Empirical Orthogonal Functions (EOFs) of the combined fields of 
near-equatorially-averaged 850 hPa zonal wind, 200 hPa zonal wind, and 
satellite-observed outgoing long-wave radiation (OLR) data. Projection of the 
daily-observed data onto such multiple-variable EOFs, with the annual cycle and 
components of inter-annual variability removed, yields principal component (PC) 
time series that vary mostly on the intra-seasonal time scale of the MJO only. 
This projection thus serves as an effective filter for the MJO without the need for 
time filtering, making the PC time series an effective index for real time use. For 
the observations, the OLR data are measured by the NOAA polar-orbiting 
satellites and is used by Wheeler and Hendon, while the NCEP/NCAR 
Reanalysis and the Global Assimilation Prediction (GASP) analyses provide wind 
fields. The index is usually available in near real time about 12 hours after the 
end of each Greenwich day (Wheeler and Hendon, 2004). 
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200—hPa Velocity Potential Anomaly: 5N-5S 



Figure 5. From mid-May into early August, weak to moderate Madden-Julian 
Oscillation (MJO) activity was observed as velocity potential anomalies 
propagated eastwards. . 
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IRMM1.RMM2) phase space for 20-Jun-2007 to 29-Jul-2007 



Labelled dots for each day. ril II I I Whsaler and Hendon (2004) 

Blue line is for Jul, green line is for Jun. BMRC Climate Forecasting 


Figure 6. The MJO index indicates that the MJO signal has weakened rapidly. 
July is indicated in blue and June in green. Courtesy of 
http://www.bom.qov.au/bmrc/clfor/cfstaff/matw/maproom/RMM/index.htm 1 

December 2007. 


2. The Local Environment of Usagi 

Deep convection in the developing system may already have started the 
transition to becoming a warm-core barotropic system even before 25 July (next 
section) due to hot towers excavating the negative vorticity and importing positive 
vorticity that then undergoes filamentation in the tail end of the baroclinic system. 

A key process is the horizontal shearing of such vorticity patches by the 
incipient vortex with an upgradient transport of vorticity into the vortex and a 
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downgradient transport into surrounding vorticity filaments. As the dynamics of 
vortex Rossby waves have been referenced to usefully characterize the wave 
mean flow interaction for the PV anomaly amplitudes considered here, it would 
also be of interest to observe wavelike features in the PV field that are consistent 
with their local dispersion relation derived in Montgomery and Kallenbach (1997). 

The filamentation, which is a consequence of the differential velocity field 
(i.e., shear, strain rate) due to one vortex acting on another that increases rapidly 
as the vortex separation distance diminishes (McWilliams, 2006). This filament 
area of high PV gradient adds to the filamentation area. The filamentation 
increases in the spiral bands, which helps to increase the circulation inside the 
intensity of the storm. The cyclogenesis of the Usagi event in the western North 
Pacific is not consistent given the scenario of TTs in the Atlantic and Gulf of 
Mexico region. In the Usagi case, aside from the broad area of anticyclone there 
was a streamer of PV anomalies pulling away from the main low depository from 
the baroclinic system. 

However, the biggest difference in the Usagi case is that the PV streamers 
detaches first from the main PV reservoir in the baroclinic system before 
undergoing cyclogenesis. The traditional paradigm for tropical cyclone formation 
in the Atlantic and southwest Pacific states “that a baroclinic disturbance in the 
mid-latitudes undergoing transformation of disturbances of sub-tropical or 
extratropical origin into tropical cyclones are often synonymous with a transition 
from a cold-core disturbance to a warm-core disturbance” (Elsberry, 2007). Can 
this occurrence happen in the western North Pacific? The Usagi transformation 
occurred while the system is already warm core as it is already below an area of 
the large anticyclone west of the decaying baroclinic system. In addition, this 
study will provide new insights into the genesis sequences of the WNP basin. 
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3. Lower Tropospheric Interactions within the Non-Developer 
(Gyre) and Developer (Usagi) 

This non-developer (gyre) initial area of convection is closely analyzed in 
the following sections and will be compared to the developer (Usagi) case. The 
850 hPa streamlines overlayed with the MTSAT imagery will highlight areas of 
convection zones of confluence. When compared with the following sections, 
you will see that the vorticity and convection difference between the two will be 
consistent with the OW and SF analysis in Section C. The investigation will show 
that the SF values in the non-developer are far less than that of the Usagi case 
as can be seen on Figure 42 and 43. The difference between the non-developer 
and the Usagi case is shown as a strong line of confluence and convection right 
along the line of confluent region in Usagi. 



Figure 7. Streamlines of horizontal flow at 850 hPa and MTSAT satellite for 0000 
UTC 22 July. Shadings are temperature values from -20°C to less than -70°C. 


At 0000 UTC 22 July, the initial non-developing system begins to detach 

the decaying baroclinic system as the anticyclone spins along behind the system 
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and stretches the convection (Figures 7). There is an area of weaker vorticity, 
moisture, and a mid-level trough sheared and damped the convection as shown 
(Figure 7). 

At 0600 UTC 22 July (not shown), the area of persistent convection is still 
being detached from the main baroclinic system. However, as the convection 
detaches a trough develops behind the anticyclone (Figure 7) This feature would 
soon damp the convection all the way thru the life cycle of this initial system. 

At 1200 UTC 22 July (not shown), the initial system completely detaches 
the baroclinic system and undergoes persistent convection. However, the large 
mid-level trough behind the developing convection is still persistent and will 
eventually shear the convection. 

By 1800 UTC 23 July (Figure 8), the initial system goes through the 
diurnal cycle and does not change in intensity as can be seen on the satellite 
animation. The system is completely detached and farther from the baroclinic 
system. As the persistent convection moves farther away from the baroclinic 
source region, the anticyclone never aligns above the convection and the mid¬ 
level trough axis clearly seen on Figure 8 will shear convection. 
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Figure 8. Streamlines of horizontal flow at 850 hPa and MTSAT satellite for 1800 
UTC 23 July. Shadings are temperature values from -20°C to less than -70°C. 


At 1200 UTC 24 July, a short wave trough (Figure 9) begins to emerge 
from the northwest and interact with the initial system (non-developer). At 1800 
UTC 24 July thru 1800 UTC 25 July, the initial system and the short wave trough 
interact and eventually enhance deep convection. The anticyclone never aligns 
to enhance upper level exhaust mechanism. There is significant increase in 
moisture as the two systems combines and move farther south into an area of 
higher SSTs, high ocean heat content, and little to no wind shear (not shown). 
This initial system is further studied in the next two sections, as to why it never 
develops. 
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Figure 9. Streamlines of horizontal flow at 850 hPa and MTSAT satellite for 1200 
UTC 24 July. Shadings are temperature values from -20°C to less than -70°C. 
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Figure 10. Streamlines of horizontal flow at 850 hPa and MTSAT satellite for 0600 
UTC 25 July. Shadings are temperature values from -20°C to less than -70°C. 


At 0600 UTC 25 July (Figure 10), the remnant of Cosme is visible to the 
east and continues to move east-northeast to be 680 nautical mile east- 
southeast of Wake Islands. Banding convection has continued to organize over 
the last 6-hours. The 850 hPa streamline at continues to indicate moderate, 
westerly vertical wind shear and this remains the primary limiting factor for 
development. Maximum sustained winds are estimated at 15 to 20 knots. 
Minimum sea-level pressure is estimated to be near 1006 hPa. With the forming 
of an area of persistent convection, the incipient (Usagi) begins to detach from 
the decaying baroclinic system. The subtropical ridge and a decaying baroclinic 
system to the north that propagates southeastward from the mid-latitudes of Asia 
dominated the synoptic pattern. This baroclinic system eventually weakened the 
subtropical ridge enough to allow the westward movement of the persistent area 
of convection (Usagi). A large anticyclonic circulation dominated the circulation 
over the eastern North Pacific over the upper mid-latitudes. 
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The animated infrared satellite depicts degraded deep convection and 
more disorganized low-level circulation center (LLCC). This weakening 
disturbance has occurred over the past 6 hours due to moderate to high amount 
of vertical wind shear as indicated by strong winds in the 850 hPa streamline and 
satellite. Although the system is moving to an area of high ocean heat content 
(not shown), the vertical wind shear is unfavorable for development at this time 
(Gray, 1968) 



Figure 11. Streamlines of horizontal flow at 850 hPa and MTSAT satellite for 1200 
UTC 25 July. Shadings are temperature values from -20°C to less than -70°C. 


At 1200 UTC 25 July (Figure 11), the tail end of the baroclinic system is 
still part of a frontal boundary and the remnant of Cosme slowly moving eastward 
roughly along 20°N as it crosses longitude 170°E. A connection between the 
outflow of Usagi and the local jet maximum associated with the mid-tropospheric 
cyclones to the northeast is not evident. The persistent area of convection 
(Usagi) continued to propagate to the southwest as the subtropical ridge 
continued to be weakened by the cyclonic circulation moving to the east. The 
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strong and large anticyclonic ridge over the mid-latitudes had advanced to be 
east of the Sea of Japan. As the persistent area of convection (Usagi) 
convection detaches from the source baroclinic system, it undergoes fluctuations 
during the diurnal cycle. This eventually limits the mechanism for deep 
convection and allows strong dry air entrainment that slows cyclone 
development. 



Figure 12. Streamlines of horizontal flow at 850 hPa and MTSAT satellite for 1800 
UTC 25 July. Shadings are temperature values from -20°C to less than -70°C. 


At 1800 UTC 25 July (Figure 12), the interaction between the outflow from 
the persistent area of convection (Usagi) and the remnants of Cosme is not 
evident. The outflow from the persistent area of convection (Usagi) has a 
pronounced anticyclonic curvature as it extends toward the mid-tropospheric low 
to the east. The pronounced anticyclonically curved outflow from the persistent 
area of convection (Usagi) associated with the mid-tropospheric low to the east 
coincides with a northward extension of the subtropical ridge. However, satellite 
reveals a diurnal variation from the previous satellite images. Tropical convective 
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systems, especially those at sea, undergo large diurnal variations because of the 
differences in the shortwave and long-wave radiation budgets of clear and cloudy 
regions (McBride and Gray, 1980). 

The diurnal variability of satellite-observed cloudiness has been described 
for the western North Pacific (Zehr, 1992). The times of maximum and minimum 
cloud area in a satellite image depend upon the cloud type and therefore the 
range of IR brightness temperatures (Tb) being considered: 

1. The strongest diurnal variation is in very cold temperature brightness 
(Tb < -65°C) clouds associated with deep convection. The maximum area of 
such cold cloud occurs between 0300-0600 Local Solar Time (LST) and the 
minimum between 1200-1800 LST. 

2. The area covered by cirrus clouds (-15°C < Tb < -45°C) also follows a 
strong diurnal cycle but with an area maximum around 1800 LST, at the time of 
the minimum in deep convective cloud area. 

3. The diurnal cycle in total cloudiness (Tb < -15°C) is weaker and follows 
that of cirrus, with a late afternoon (1800 LST) maximum. 

4. Heavy rain follows the same diurnal cycle as convection. The 
maximum occurs between local midnight and dawn over the ocean, but can be 
anytime over land, especially when orographic effects dominate. 

Tropical cyclone structure changes, on the other hand, show little 
substantiated diurnal preference. Because many forecasting methods use trends 
in satellite imagery, it is critical that forecasters be aware that cloudiness, and 
therefore the satellite signatures of tropical systems, also show a marked diurnal 
variation unrelated to structural changes. Short-term trends in satellite imagery 
can be very misleading and cause major forecast errors (BOM 2008). 
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Figure 13. Streamlines of horizontal flow at 850 hPa and MTSAT satellite for 0000 
UTC 26 July. Shadings are temperature values from -20°C to less than -70°C. 


At 0000 UTC 26 July (Figure 13), the subtropical ridge weakened enough 
to allow Usagi to begin recurvature to the southwest. The trough (old baroclinic 
system) to the northeast of Usagi became nearly fully aligned with the tropical 
cyclone. The continued southwestward motion of Usagi initiates the merger 
between the outflow from the tropical cyclone and the mid-latitude anticyclone to 
the north. The outflow region of the persistent area of convection (Usagi) 
continues to exhibit pronounced anticyclonic curvature in association with the 
mid-tropospheric baroclinic system to the northeast. 
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The outflow of the persistent area of convection (Usagi) and the 
anticyclone to the northwest have combined to form a single wind maximum 
exceeding 50 ms' 1 north of the tropical cyclone. However, the downstream 
outflow continues to curve anticyclonically toward the mid-tropospheric cyclone, 
while the mid-latitude flow maintains a meridional orientation. Furthermore, the 
subtropical ridge continues to build northward between the persistent area of 
convection (Usagi) and the 850-hPa low associated with baroclinic system to the 
east. 



Figure 14. Streamlines of horizontal flow at 850 hPa and MTSAT satellite for 0600 
UTC 26 July. Shadings are temperature values from -20°C to less than -70°C. 
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At 0600 UTC 26 July (Figure 14), the persistent area of convection (Usagi) 
and Cosme have become close enough that if there would be interaction this 
would be the time. However, as analysis suggest they are still at least 10° apart 
and too far to form interaction as mentioned earlier. The northward advection of 
warm air is defined as the anomalously warm regions (not shown) to the west of 
the storm center. At low levels, warm anomalies to the east (not shown) are 
associated with the warm-air advection. A mid-level warm anomaly is found over 
the center in association with the developing tropical cyclone (not shown). 



Figure 15. Streamlines of horizontal flow at 850 hPa and MTSAT satellite for 1200 
UTC 26 July. Shadings are temperature values from -20°C to less than -70°C. 


At 1200 UTC 26 July (Figure 15), the persistent area of convection (Usagi) 
continued to recurve and became fully imbedded in the anticylonic ridge 
propagating off the mid-latitudes. The relative positioning of the tropical and mid¬ 
latitude circulations was indicative of the start of the formation stage of SI. 


40 









Figure 16. Streamlines of horizontal flow at 850 hPa and MTSAT satellite for 0000 
UTC 27 July. Shadings are temperature values from -20°C to less than -70°C. 


At 0000 UTC 27 July (Figure 16), the persistent area of convection (Usagi) 
outflow and anticyclone to the northwest have closely merged. The decaying 
baroclinic system is now fully washed out to the east in conjunction with the 
forcing of vertical motion north of the persistent area of convection (Usagi). This 
vortex interaction as seen on VORTRACK (Figure 2-4) could have started about 
this time. However, the weaker vortex that is the remnants of Cosme being 
advected poleward from the stronger baroclinic vortex, which is the tail end of the 
baroclinic system. 


41 







Figure 17. Streamlines of horizontal flow at 850 hPa and MTSAT satellite for 0600 
UTC 27 July. Shadings are temperature values from -20°C to less than -70°C. 


At 0600 UTC 27 July (Figure 17), the area of persistent deep convection 
(Usagi) is becoming consistent near 19°N 155°E. This is approximately 570 
nautical miles east northeast of Saipan. The animated infrared imagery indicates 
a rapidly consolidating LLCC with deep convection. The disturbance is 
imbedded in a tail end of a baroclinic system over the central Pacific. The 850 
streamline analysis suggest a large anticyclone west of the decaying baroclinic 
system that will enhance equator ward flow and low vertical wind shear ideal for 
deep convection and will move this area of deep convection to a favorable 
environment for cyclogenesis. This rapid consolidation of LLCC and favorable 
region for convection enhances this tropical seedling to develop and form a 
significant tropical cyclone. 
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Figure 18. Streamlines of horizontal flow at 850 hPa and MTSAT satellite for 1200 
UTC 27 July. Shadings are temperature values from -20°C to less than -70°C. 


At 1200 UTC 27 July (Figure 18), the deep convection results in a strong 
cyclonic curvature toward the mid-tropospheric anticyclone to the east of Usagi. 
The cyclonic curvature is maintained as the persistent area of convection (Usagi) 
moves northward and the mid-tropospheric cyclone to the east digs southward at 
1800 UTC 27 July (Figure 19). 
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Figure 19. Streamlines of horizontal flow at 850 hPa and MTSAT satellite for 1800 
UTC 27 July. Shadings are temperature values from -20°C to less than -70°C. 
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Figure 20. Streamlines of horizontal flow at 850 hPa and MTSAT satellite for 0600 
UTC 28 July. Shadings are temperature values from -20°C to less than -70°C. 


At 0600 UTC on 28 July (Figure 20), the SI has begun by increased 
moisture and increased cyclonic circulation of the system. The system also has 
moved farther south into an area of higher SSTs, high ocean heat content, and 
little to no wind shear. Thus, the remnants of the tail end of the baroclinic system 
eventually became Typhoon Usagi. Usagi initially moved southward in response 
to an anticyclone to its northwest, and then northwestward toward Japan after 
having as it moved around the southern periphery of the ridge. 
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Figure 21. Streamlines of horizontal flow at 850 hPa and MTSAT satellite for 1200 
UTC 28 July. Shadings are temperature values from -20°C to less than -70°C. 


At 1200 UTC on 28 July (Figure 21), as TD 5W (Usagi) entered the 
formation stage of SI a clear interaction occurred between the outflow of the 
decaying baroclinic wave circulation and a mid-tropospheric anticyclone to the 
west. The outflow first curved anticyclonically toward the cyclone to the east. 
The curvature became pronounced as Usagi moved northward and the cyclonic 
circulation drifted southward. This interaction coincided with a pronounced 
northward extension of the subtropical ridge. In this case, TD 5W (Usagi) was 
classified as a tropical depression 05W at 1200 UTC 28 July, which defines 36-h 
duration for the transformation stage (JTWC, 2007). 
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At 0000 UTC on 29 July (not shown), TS Usagi (5W) entered its strongest 
stage to date as an area of deep convection where the upper anticyclone is 
aligned for the much-needed exhaust of deep convection. At this time of SI, a 
clear interaction occurred between the outflow of the decaying baroclinic wave 
circulation and a mid-tropospheric anticyclone to the west. 

At 0000 UTC on 30 July (not shown), Usagi was named as Typhoon Usagi 
(5W) and is now 315 nautical miles south-southeast of Iwo-Jima, tracked 
northwestward at 11 knots. Maximum sustained winds were estimated at 65kt 
gusting to 80kt (JTWC, 2007). 

B. DYNAMIC AND THERMODYNAMIC EVOLUTION OF USAGI 


This section will investigate high-resolution analysis of the GFS FNL 
model to show that the PV of Cosme is too far away from the baroclinic system 
and did not play any essential role in the development of Usagi. The analysis will 
show how Cosme could have aided the development of Usagi via an indirect 
factor. As mentioned in the previous section, the apparent vortex merger starts 
about 0000 UTC on 27 July as shown in figure 2, and with each passing chart, 
the weaker vortex that is the remnants of Cosme are being advected to the 
stronger baroclinic vortex, which is the tail end of the baroclinic system and 
eventually became Usagi. The PV from the baroclinic disturbance could have 
become close enough to wrap around the filaments from Cosme into the 
disturbance that eventually became Usagi. 

Potential vorticity is defined on constant potential temperature isobaric 
surfaces in layer average between 500 and 850 hPa (Holton, 2004) as follows: 


PV = -g 


dV 


86 _ 

dp 


where 


( 2 . 1 ) 


g is the gravitational acceleration, 
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£ is the relative vorticity on a pressure surface, 


f is the Coriolis parameter, 

V is the layer-average vector wind at 500 and 850 hPa, 
0 is potential temperature, 


and p is pressure. 

The PV anomaly is defined by the deviation from the zonal average PV. 
For adiabatic motions potential vorticity is a materially conserved quantity. It is 
not materially conserved when moist processes are active. Using synoptic-scale 
gridded fields, the computation of the diabatic influences can be difficult to 
compute. However, the total derivative of PV can be used to define a proxy for 
diabatic processes as: 


d(PV) 

dt 


P 


where 


( 2 . 2 ) 


p is the density, 


@ is the diabatic heating, 

and the overbar represents the layer average. 

A tropical cyclone is represented by a large deep vertically upright region of 
positive PV. As the seeding PV anomaly undergoes SI, warm advection begins 
in the eastern quadrants of the storm, while cold advection occurs in the western 
quadrants. The warm, moist air that is advected around the western side of the 
storm is associated with rising motion. The warm advection and upward vertical 
motion and associated release of latent heat will contribute to an increase in 
thickness to the east of the southward-moving tropical cyclone. In terms of the 
distribution of PV, the mid-tropospheric heating maximum on the western side of 
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the decaying baroclinic cyclone will contribute to a positive change in the lower- 

d{Pv) d {Pv) 

level PV ( dt >0), and a negative change in the upper-level PV ( dt <0). 
At the same time, the warm advection occurring on the eastern side of the storm 
brings warm, moist air into the storm, which enhances the convection and 
eventually the eye wall convection. This enhancement eventually causes the 
storm to take on the more symmetric appearance of a tropical cyclone. To 
diagnose the dynamic evolution of the PV distribution, plan-view diagrams of the 
PV, PV anomaly, the time derivative of the potential vorticity, relative humidity 
and vertical motion (u>) are constructed from the gridded analyzed fields. 

The time tendency of the PV distribution is used to diagnose areas where 
diabatic heating is influencing the PV distribution at upper and lower levels 
(Vancas 2007). 

Relative humidity (SF) and vertical motion (OW) are used to indicate the 
potential for sustained deep convection. 

Tropical cyclone formation brings regions of anomalously warm 
temperatures, particularly in the lower atmosphere, poleward. For a baroclinic 
system undergoing decay, warm anomalies typically have an upright structure 
similar to that of the PV distribution discussed above. The peak warm anomaly 
in this case occurs in the lower atmosphere. The amount of the anomalous 
warming may be related to the intensity of the tropical cyclone, with stronger 
storms having higher peak anomaly temperatures (Krishnamurti et. al. 1987). 
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Figure 22. Potential vorticity in units of PVU at 850-500 hPa overlay for 1800 UTC 
25 July 


By 1800 UTC 25 July (Figure 22), the distribution of PV is becoming 
evident as a the persistent area of convection (Usagi) and the approaching 
remnants of Cosme. However, regions of positive and negative PV tendency 
(not shown) are beginning to appear centered near 850-500 hPa to the east and 
west of Usagi. 
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Figure 23. Potential vorticity in units of PVU at 850-500 hPa overlay for 0000 UTC 
26 July 


By 0000 UTC 26 July (Figure 23), the region of positive PV had increased 
in the mid-troposphere on the persistent area of convection (Usagi). The 
distribution of OW values becomes stronger than it was 6-hour earlier (Appendix 
A). 
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Figure 24. Potential vorticity in units of PVU at 850-500 hPa overlay for 0600 UTC 
26 July 


By at 0600 UTC 26 July (Figure 24), the PV has become concentrated 
over the storm center at low levels. The vertical distribution of upward vertical 
motion has changed dramatically. Sinking motion is found west and east near 
the center of Usagi. Therefore, between 1800 UTC 25 July and 0000 UTC 26 
July the circulation of Usagi seems to have become better defined. This may be 
in response to a period of reduced vertical wind shear. Strong vorticity at the tip 
of the PV streamer favors deep convection, which will enhance the low-level 
vorticity and moisten the mid level, while the remnants of the Cosme wave may 
act to enlarge the wave “pouch” and help protect the moisture in the middle 
levels from outside influences. In this sense, the Cosme wave may have played 
a role in the development of the persistent area of convection (Usagi) even 
though it did not directly interact with the baroclinic system and Usagi. At this 
time, the region between Usagi and Cosme also reveals a small area of clearing 

between the two consistent areas of convection. 
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Figure 25. Potential vorticity in units of PVU at 850-500 hPa overlay for 1200 UTC 
26 July 


By 1200 UTC 26 July (Figure 25), several important features can be 
noted. It is from this time in the SI event that a more detailed description will be 
given in the following sections when we review the OW and SF analysis. The 
anticyclone is evident behind the baroclinic system to the west of Usagi (as seen 
in the previous section) as a lobe of low PV and an area clearing behind it. 
Furthermore, an area of positive at low levels just to the west of the storm center 
is associated with the anticyclone that is a broad region of subsiding air to the 
west of the center of Usagi. Because the lower troposphere warm anomaly (not 
shown) remains over the storm center, there is also a southern extension of low 
level warming. In the lower levels, the vorticity anomaly becomes organized with 
enhanced low-level convergence from the existing instability from the low-level. 
This would be consistent with the bottom-up formation of a frontal system 
stagnating in the subtropics and transforming to warm-core characteristics (DMW 
2008). 
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Figure 26. Potential vorticity in units of PVU at 850-500 hPa overlay for 1800 UTC 
26 July 


By 1800 UTC 26 July (Figure 26), the area of low-level positive PV to the 
center Usagi extends to the mid-troposphere as the mid-latitude anticyclone 
approaches the forming storm (Usagi). A large area of dry air (not shown) exists 
to the west of the storm center with some subsidence. However, the warm 
anomaly through the entire troposphere has been shifted to the east of the center 
of Usagi. Significant changes in the distribution of PV, time tendency of PV, and 
temperature distribution now occur in the 6 hour between 1800 UTC 26 July and 
0000 UTC 27 July. As the large anticyclone approaches Usagi from the west 
seen the previous section, the region of positive PV time tendency (not shown) is 
nearly adjacent to Usagi, and the distribution of PV in the mid-troposphere has 
begun to extend into Usagi. The area of stronger PV of Cosme moves to cooler 
waters and eventually dissipates. 


54 






Figure 27. Potential vorticity in units of PVU at 850-500 hPa overlay for 0000 UTC 
27 July 


By 0000 UTC 27 July (Figure 27), at low-levels the distribution of vertical 
motion is distinctly over to the center of Usagi. Cosme has decreased in the 
strength at this time, as its PV gets weaker and smaller in horizontal scale, the 
low-level maximum PV anomalies of Usagi are shifted southerly. 
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Figure 28. Potential vorticity in units of PVU at 850-500 hPa overlay for 0600 UTC 
27 July 


By 0600 through 1200 UTC 27 July (Figure 28 and 29), the area of 
persistent convection (Usagi) has taken on characteristics of the formation stage 
of SI. The maximum positive PV region in the low-levels has moved over the 
center of the storm as the developing circulation of the convection (Usagi) moves 
to an area of higher SSTs and weak vertical wind shear favorable for TC 
formation. 
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Figure 29. Potential vorticity in units of PVU at 850-500 hPa overlay for 1200 UTC 
27 July 


As Usagi entered the low latitudes, it underwent several important 
dynamic and thermodynamic changes. The distribution of upward vertical motion 
underwent dramatic changes under the influence of the mid-latitude circulations. 
As the SI stage was completed, a strong region of upward vertical motion 
developed farther to the east associated with the remnants of the cyclonic 
circulation of the lower tropospheric baroclinic system. A secondary upward 
vertical motion region centered east and over the decaying baroclinic system was 
associated with the system moving off the mid-latitudes. 

The temperature distribution also underwent complex changes as seen in 

the MTSAT data. As expected for a developing tropical cyclone with such 

baroclinic origin, the temperature distribution initially had the maximum warm 

anomaly in the lower troposphere centered over the storm. Under the influence 

of the mid-latitude circulations, the low-level warm anomaly first shifted to the 

east and north of the storm center. To examine the processes associated with 
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the southward extension of the baroclinic system, the detachment of the PV 
anomaly and the lower-level wind field is now examined. 

The original hypothesis that Cosme was essential in the development of 
Usagi and this study was at first motivated by this hypothesis, the analysis shows 
that there is no vorticity connection between the two systems, there an area of 
subsidence between the baroclinic system and the remnants of Cosme on 0600 
UTC 26 July when both at are close proximity to each other (10° apart). If indeed 
Cosme aided the formation of Usagi it may have been through another means 
rather than direct vorticity transfer between ex-hurricane Cosme and pre-Usagi. 
The analysis suggests that the Cosme wave may have played an important 
indirect role in the formation of Usagi by enlarging the wave pouch and helping 
protecting the proto-vortex from hostile outside influences. Further evaluation of 
this theory will be shown in the next section by comparing OW, relative humidity, 
and SF fields. This is confirmed using VIS-5D simulations (Figures 30-32). The 
VIS-5D illustrated that there is no vorticity or trajectory connection between the 
two systems. One method to further explore the direct influence of ex-Cosme 
would be to remove the Cosme circulation from the initial condition of a model 
and compare the predictions with and without Cosme circulation; however such 
an experiment is beyond the scope of this study. 
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Figure 30. Vis5D trajectories between Usagi and Cosme during 0600 UTC 26 
July. Streamlines are 850 hPa. The orange figures are absolute vorticity 
(lO^S 1 ). The white strips are trajectories with Usagi and the yellow are 
trajectories associated with Cosme from the last 72 hours. 
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Figure 31. Vis5D trajectories between Usagi and Cosme during 1200 UTC 27 
July. Streamlines are 850 hPa. The orange figures are absolute vorticity 
(lO^S 1 )). The white strips are trajectories with Usagi and the yellow are 
trajectories associated with Cosme from the last 72 hours. 


60 





Figure 32. Vis5D trajectories of Usagi during 1800 UTC 28 July. Streamlines are 
850 hPa. The orange figures are absolute vorticity (lO^S 1 ). The white strips are 
trajectories with Usagi and the yellow are trajectories associated with Cosme 
from the last 72 hours. 

C. OW CRITERION AND SF PARAMETERS TO DISTINGUISH THE NON 
DEVELOPER (GYRE) AND DEVELOPER (USAGI) 

In Chapter I, the relative vorticity and OW criterion are discussed as they 
are used to distinguish curvature vorticity from shear vorticity. The large positive 
value of OW parameter suggests regions where curvature vorticity is dominant. 
The moving frame of reference is optimal to display the Lagrangian flow and the 
evolution of the storm vortex. Isolating these parameters, you can see true 
Lagrangian parcel motions that are affecting the dynamics and thermodynamics 
of the transformation process in the ground based frame of reference. Here, we 
use the criterion based on the OW parameter to extract the vortex statistics from 
the full vorticity field. 
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This section will test the hypothesis that the tip of the PV streamer, which 
is dominated by curvature vorticity, is the key factor leading to the TC formation. 
Strong vorticity may be sufficient to initiate deep convection, which will then in 
turn enhance the low-level vorticity and moisten the middle level. 

The distributions of vertical motion (PV) and moist air (SF) suggest that 
the relationship of the thermodynamics along the tip of the decaying baroclinic 
system contributes to a mid-troposphere moist anomaly that will be confirmed by 
TRMM data. The relative humidity SF parameter steadily increases (in percent) 
beginning at 1200 UTC 25 July and reaches its strongest value (most moist) at 
0000 UTC 27 July. 

The dynamic and thermodynamic processes discussed in Chapter III.B 
are compared with the relative vorticity, OW, and SF findings, respectively, and 
appear to be important precursors associated with the SI the non-developer and 
Usagi. The dynamic and thermodynamic forcing of these parameters evidence 
noticeable changes in the middle and lower-troposphere are noticeable in 
changes in OW, SF and PV anomaly magnitude. 

1. OW Parameter Comparison 

By comparing the two disturbances, using previous finding from PV, 
relative vorticity and OW values, the following interesting features distinguished 
the non-developer from the developer. 

The PV streamer within the baroclinic system has a positive value of OW 
at its tip, which suggests that vorticity is dominant in that area. This became 
detached from the PV streamer and subsequently formed the non-developer and 
Usagi (Figures 33 and 34). At 1800 UTC 23 July, the initial system detaches 
from the decaying baroclinic system as anticyclone spins along behind the 
system and stretches the convection. This persistent convection starts to interact 
with an existing trough axis behind the baroclinic system. This area of 
convection will become persistent as it moves thru a region favorable for TC 
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formation. However, as the persistent area of convection detaches from the 
main baroclinic region, a mid to upper-level trough behind the anticyclone begins 
to develop (Figure 33). This feature would soon shear and damp the convection 
all the way thru the life cycle of this initial system. Unfortunately, as the 
persistent convection of the initial system moves farther away from the baroclinic 
source region, the anticyclone never aligns above the initial system as the 
impeding trough axis damps the cyclogenesis process. Thus, this initial system 
never develops. There is significant increase in moisture as the two system 
combines and move farther south into an area of higher SSTs (not shown), high 
ocean heat content (not shown), and little to no wind shear (not shown). 
However, with the necessary and ideal component for deep convection, this 
initial system never develops into a tropical cyclone. 

In contrast with the non-developing system, the Usagi case exhibits a 
different scenario. As the Usagi persistent area of convection detaches the main 
baroclinic source region, the PV streamers within the baroclinic system is as 
strong as the initial system and has a positive value of OW at its tip, which 
suggests that vorticity again is dominant in that area as can be seen in the four- 
panel evolution of Usagi (Figure 34). As it moves towards southeast, an area 
favorable for deep convection, the baroclinic system is already stronger near 
surface and the approaching remnants of wave Cosme is stronger at 600 and 
850 hPa. During 1200 UTC 27 July, at 600 and 850 hPa (Figures 34 and 36) in 
the frame of reference moving with the baroclinic system, a big pouch formed as 
Cosme approached the baroclinic system and two sub-gyres (Usagi and Cosme) 
formed within the mother pouch. The one in the west (Usagi), which is directly 
related to the baroclinic system later became TD 5W on 1200 UTC 28 July. 
Although the baroclinic system is stronger near the surface at 850 hPa, there 
was no closed circulation in the moving frame of reference until 6 hours before 
genesis at 1200 UTC 27 July as can be seen at 600 hPa. 

By comparing the non-developer (gyre) and developer (Usagi) based on 
the four panel 600 and 850 hPa OW and relative vorticity analysis (Figure 33-40), 
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the non-developer has almost the same values of vorticity compared to the 
developer. However, the non-developer encountered environmental limitations 
such as a mid to upper-level trough axis interaction and lack of protection from 
outside elements that eventually limit the internal mechanisms for deep 
convection, allowing for strong dry air entrainment that inhibits cyclone 
development.. These processes inhibits deep convection needed for TC 
formation, thus the initial invest region never develops. Whereas, the Usagi 
system had environmental protection in the critical level needed for deep 
convection persistent enough for TC formation as it moves to a favorable region 
of cyclogenesis. As introduced in Chapter I, the marsupial paradigm discusses 
the proto-vortex cyclonic eddies instrumental in TC formation are intimately 
associated with parent wave’s critical latitude in the lower-troposphere. This 
structure forms a critical layer because of the wave’s finite amplitude interaction 
with its own critical latitude, is a region of cyclonic rotation and weak 
strain/shearing deformation in which synoptic and mesoscale anomalies move 
together and amplify on a nearly zero relative mean flow. This multi-scale 
interaction provides dynamically pathway to bottom-up development of the proto 
vortex from below. The critical layer of the parent wave provides a set of closed 
material contours inside of which air is repeatedly moistened by convection. The 
wave is protected to some degree from lateral intrusion of dry air impinging 
vertical shear, and able to keep in pace with the parent wave until the proto¬ 
vortex has strengthened into a self-maintaining entity and begins control its own 
destiny. By correlating the marsupial analogy for TC formation in this study, we 
can verify in theory that the remnants of Cosme acted as a protective proto 
vortex associated with the parent’s wave (Usagi) critical latitude at 600 hPa as 
can be clearly seen on Figure 34 panel 2 as a closed streamline of protected 
area. 

The findings above would be consistent with the hypothesis that Cosme 
may not be essential in the development of Usagi, but have just probably 
widened the area of the “pouch” of vorticity anomalies. This anomaly can be 
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seen at 600 hPa on 1200 UTC 27 July (Figure 34 panel 2). This area of strong 
vorticity at the tip of the PV streamer favors deep convection, which will enhance 
the low-level vorticity and moisten the mid level, while the Cosme wave enlarges 
the wave pouch and helps to preserve moisture in the middle level of Usagi 
convection. In this sense, the Cosme wave played an important role although it 
did not directly interact with the baroclinic system. 

The relative vorticity Hovmoller diagram as seen on Figure 41 suggest a 
consistent vorticity motion of -7 ms' 1 from 1000 thru 500 hPa. This is consistent 
that vorticity is all the way down in the low-levels. The Usagi storm propagation 
is nearly parallel to the trough propagation. 
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Figure 33. Streamlines of horizontal flow at 600 hPa for the genesis sequence of 
the non-developer from 1800 UTC 23 July to 1200 UTC 24 July. The shadings 
indicate high values of OW parameter (units: 10' 10 s' 2 ) as defined in Eq. (1.1). 

This quantity, like vorticity, is invariant with respect to translation, therefore 
identical in translating and moving frames. Red areas are strong values and blue 
values areas are weak. 
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Figure 34. Streamlines of horizontal flow at 600 hPa for the genesis sequence of 
the Usagi system from 0600 UTC 27 July to 0000 UTC 28 July. The shadings 
indicate high values of OW parameter (units: 10' 10 s' 2 ) as defined in Eq. (1.1). 

This quantity, like vorticity, is invariant with respect to translation, therefore 
identical in translating and moving frames. Red areas are strong values and blue 
values areas are weak. 
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Figure 35. Streamlines of horizontal flow at 850 hPa for the genesis sequence of 
the non-developer from 1800 UTC 23 July to 1200 UTC 24 July. The shadings 
indicate high values of OW parameter (units: 10' 10 s' 2 ) as defined in Eq. (1.1). 

This quantity, like vorticity, is invariant with respect to translation, therefore 
identical in translating and moving frames. Red areas are strong values and blue 
values areas are weak. 
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Figure 36. Streamlines of horizontal flow at 850 hPa for the genesis sequence of 
the Usagi system from 0600 UTC 27 July to 0000 UTC 28 July. The shadings 
indicate high values of OW parameter (units: 10' 10 s' 2 ) as defined in Eq. (1.1). 

This quantity, like vorticity, is invariant with respect to translation, therefore 
identical in translating and moving frames. Red areas are strong values and blue 
values areas are weak. 
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Figure 37. Streamlines of horizontal (rotational + divergent) flow and Zeta at 600 
hPa of the non-developer from 1800 UTC 23 July to 1200 UTC 24 July. Shading 
indicates relative vorticity (units: 10' 5 s' 2 ). 
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Figure 38. Streamlines of horizontal (rotational + divergent) flow and Zeta at 600 
hPa of the Usagi system from 0600 UTC 27 July to 0000 UTC 27 July. Shading 
indicates relative vorticity (units: 10' 5 s" 2 ). 
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Figure 39. Streamlines of horizontal (rotational+divergent) flow and Zeta at 850 
hPa of the non-developer from 1800 UTC 23 July to 1200 UTC 24 July. Shading 
indicates relative vorticity (units: 10' 5 s' 2 ). 
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Figure 40. Streamlines of horizontal (rotational+divergent) flow and Zeta at 850 
hPa of the Usagi System from 0600 UTC 27 July to 0000 UTC 28 July. Shading 
indicates relative vorticity (units: 10' 5 s" 2 ). 
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Figure 41. Relative vorticity Hovmoller diagram 3 days before development of 
Usagi from 26 -30 July. 
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2. SF Comparison 


The time evolution of SF and RH at 600 hPa Figure 42 and 43 in the 
moving frame of reference of the baroclinic system) was examined and the two 
disturbances, the non-developer and developer were compared. As mentioned 
in Chapter 1, the SF of the lower troposphere, defined as the ratio of precipitable 
water to saturated precipitable water. 

By comparing the two disturbances using SF values, the non-developer 
over the developer, it is suggested that the non-developer had less precipitable 
water (lower SF) than the developer (Usagi). As discussed in the previews 
section, the non-developer also became detached from the source baroclinic 
system, it undergoes fluctuations as it gets under the influence of an mid-level 
trough (can be seen on Figure 42 as an area of lower SF values west of the non¬ 
developing system) as it moves southeastward thru the lower latitudes. The non¬ 
developer is moving to an area of low vertical wind shear that is ideal for deep 
convection and moves into an area that is favorable environment for 
cyclogenesis area of deep convection. In comparison, using the precipitation 
distribution it is evident that the TRMM distribution (Figures 44 and 45) is 
consistent with SF analysis. The low SF values eventually limit the mechanism 
for deep convection, allows strong dry air entrainment that slows cyclone 
development to weaken, and then the system becomes disorganized (Elsberry 
2007). 

In contrast with the non-developer, major changes in the thermodynamic 
structures of the developer (Usagi) began at 1200 UTC 26 July due to the 
influence of the low-level anticyclone to the west and the high distributions of 
temperature (SST), moisture (SF), and large values of OW. The maximum OW 
anomaly and SF parameters in the 600-hPa is consistent with the TRMM that 
occurs in a round-shaped region over the northern half of the circulation with a 
maximum to the west of the storm center (SF compared to TRMM four panel). 
The upward vertical motion seen as vorticity (OW) to the storm center coincides 
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with a broad, low-level PV anomaly seen in the entire PV figure in the past 
section. A relative moist anomaly occurs to the storm center at 600 hPa is 
consistent with the SF values (Figure 42 and 43). The SF evolution and 
translating streamlines in the development of Usagi is shown in Figure 43 at 600 
hPa. The closed streamlines of directional flow are formed near the trough axis, 
delineating a recirculating gyre of finite width. Strong cyclonic/rotational vorticity, 
high SF and weak vertical wind shear simultaneously coincide within the gyre 
and provide a favorable sub-synoptic scale environment for TS formation. The 
vertical signature of the Cosme wave and diabatic vortex prior to Usagi formation 
is clearly consistent with TRMM data suggesting that convective activity is 
present near the critical latitude 24 hours prior official TD 5W (Usagi) formation in 
1200 UTC 28 July (JTWC 2007). Near the critical latitude, the convective activity 
near pre-Usagi becomes sustained and more focused, while the Cosme wave 
moves northeastward away from the critical latitude and dies out, consistent with 
the conclusions of DMW08. 

Compared to the non-developing intrusion, the stronger Usagi has 
consistently high values of SF as it detaches from the baroclinic system, 
approaches the remnants of Cosme, and protrudes into the lower latitudes. The 
pre-Usagi system is stronger near the surface. The Cosme wave has become 
less moist (<50%) as it gets near with Usagi by 0600 UTC 27 July (Figure 43 
Panel 2). However, at 600 hPa in the frame of reference moving with the 
baroclinic system, a pouch formed as Cosme approached the baroclinic system 
and two sub-gyres formed within the parent pouch. This area of high SF values 
favors deep convection, which will enhance the low-level vorticity and moisten 
the mid level, while the Cosme wave appears to enlarge the wave pouch and 
helps moisten the middle levels. In this sense, again a hypothesis, Cosme wave 
played an important role although it did not directly interact with the baroclinic 
system. The convective activity that is Usagi becomes sustained and more 
focused, while the Cosme wave moves northeastward away from the critical 
latitude and dies out. 
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The Usagi system formed in relatively moist, precipitating environment 
with large values of vorticity as seen with the PV and OW parameter. This is true 
for a gyre and local averages alike, but the local average contains some 
exceptionally large values of precipitation and OW parameter, suggesting that 
the TRMM satellite has observed enhanced precipitation associated with the 
developing storm, concentrated at the best track location rather than distributed 
around the gyre. The ECMWF analysis have generally pinpointed the location of 
the developing vortex within the gyre (Figure 42 Panel 2) and identified its 
vorticity as predominantly rotational rather deformational, i.e. not deformed by 
strain or shear. The persistence of precipitation defined in TRMM imagery 
preceding formation is evident. The formation of closed streamlines implies that 
air with the closed space is protected, to some degree, from lateral entrainment 
of dry air. This containment is important in allowing a convective type of heating 
profile to dominate the closed region with relatively smaller contributions from 
stratiform precipitation that are otherwise common in mesoscale convective 
systems (Tory and Montgomery, 2006). 
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Figure 42. Streamlines of horizontal flow 600 hPa for the genesis sequence of the 
non-developer, with high values of saturation fraction from 1800 UTC 23 July to 
1200 UTC 24 July. Red areas are strong values and white areas are weak 
values of SF. 
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Figure 43. Streamlines of horizontal flow 600 hPa for the genesis sequence of the 
Usagi system, with high values of saturation fraction from 0600 UTC 27 July to 
0000 UTC 28 July. Red areas are strong values and white areas are weak 
values of SF. 



Figure 44. Streamlines of horizontal flow 850 hPa for the genesis sequence of the 
non-developer, with high values of TRMM 3-hour accumulated precipitation 
indicated by shading from 1800 UTC 23 July to 1200 UTC 24 July. Red areas 
are strong values and white areas are weak values of precipitation. 
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Figure 45. Streamlines of horizontal flow 850 hPa for the genesis sequence of the 
Usagi system, with high values of TRMM 3-hour accumulated precipitation 
indicated by shading from 0600 UTC 27 July to 0000 UTC 28 July. Red areas 
are strong values and white areas are weak values of precipitation. 

3. Time Series Analysis 

Time series analysis was constructed to show the relationship between 
the two systems. This experiment will be used to summarize the relationships 
between saturation fraction distribution area averages. The RH values at 700 
hPa were averaged at each synoptic time from 0000 UTC 22 July until 1200 UTC 
25 July in a 10 deg. lat. by 10 deg region centered eight deg. lat. of the storm for 
the non-developer and from 0000 UTC 26 July until 0000 UTC 29 July for Usagi. 
The selection of this region allowed the detection of the parameters as the storm 
progressed southwestward. By using GARP MTSAT overlay and Vis5D absolute 
vorticity parameters, the convection calculation is as follows:. 

Average RH = (RH max + RH min)/2 

The saturation fraction (RH) was averaged over a moving frame between 15°- 
30°N and from 140°-165°E, which was the region directly on track of Usagi and 
the non-developer. For the Cosme wave, the SF is almost less that 60% as it 
interacts with the baroclinic system on 1200 UTC 27 Jul (not shown). The 
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Usagi case on the other hand has constant high values of SF at more than 85% 
all throughout formation time compared to Cosme and the non-developer (Figure 
46 and 47). 

The convective available potential energy (CAPE) time-series values 
associated with the non-developer and Usagi were calculated over the same time 
and moving frame of reference with the RH, which was the region directly on 
track of the persistent area of convection. By comparison, the non-developer has 
increasing values up to 2300 J/Kg then suddenly decreased while the Usagi case 
had increasing values up to 2700 j/Kg (Figure 48-49). A more complete estimate 
of buoyant energy can be computed by determining the temperature difference 
between the ascent path and the environment at all levels from the level of free 
convection up through the equilibrium level. This area provides a measure of the 
integrated effects of the potential temperature differences between the rising 
parcel and its environment. The quantitative measure of this positive area is 
CAPE. CAPE is measured in units of joules per kilogram (J/kg). This layer 
represents the average heat and moisture conditions fueling deep convection. 
The environmental CAPE ideal for deep convection and modulate convection is 
often in the range of 1000-2000 J/kg (METED, 2007). 

Detailed time lines and figures or OW, Zeta, TRMM, and SF of the non¬ 
developer (Gyre) and developer (Usagi) are located in the Appendix. 
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Relative Humidity Time Series of Non-Developer 



Figure 46. Relative humidity time series diagram of non-developer from 22 -25 
July. 
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Figure 47. Relative humidity time series diagram of Usagi from 26 -28 July. 
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IV. SUMMARY AND RECOMMENDATION 


A. SUMMARY 

During July 2007, two distinct areas of disturbed weather with Sis 
emerged out of remnants of a decaying baroclinic system in the WNP. The initial 
invest system became detached but never developed and the second identified 
system eventually became Typhoon Usagi (5W). The decaying baroclinic 
system with a trailing cold front pushed through the WNP and the associated 
cyclogenesis occurred along the trailing edge of the cold front. (This was 
illustrated with the areas of concentrated 850 hPa PV and relative vorticity 
moving down the front.) The second persistent area of convection (Usagi) was 
already near the surface and was rather strong, so it did not take much to convert 
this cold core system into a tropical system given the fact it was moving south 
into warm SSTs, low vertical wind shear.. 

The keys to events that led to the non-development of the initial invest 
area and the developments of Usagi were investigated. Fields from the NCEP 
GFS final analyses and ECMWF model analysis of OW and SF each in 6-hour 
intervals suggested the development of two distinct areas of disturbed weather 
with sub-tropical cyclogenesis. The Vis5D was a valuable tool for visualizing 
both gridded data and irregularly located data. Sources for these data included 
numerical weather prediction models output, surface observations, and other 
similar sources. The high-resolution model analysis was used and we concluded 
that the events leading to the Usagi TC formation is not vortex merger and 
axisymmetrization of the ex-hurricane Cosme, as initially hypothesized. 

By comparing the two sub-tropical intrusions, using regions where 

curvature vorticity is dominant (OW) and the ratio of precipitable water to 

saturated precipitable water (SF) is significantly above environmental values, it is 

found that the non-developer (gyre) has significantly less concentration and 

organization of OW and SF than the developing disturbance (Usagi). The non- 
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developer become detached from the source baroclinic system, and undergoes 
fluctuations as it interacts with a mid to upper-level pressure trough (an area of 
low SF values (< 50%) west of the non-developing system. In comparison to the 
Usagi system, the non-developer has less precipitable water and moves to lower 
latitudes due to environmental factors as it get detached and develops away from 
the source baroclinic system. The disorganized area of OW and low SF values 
are hypothesized to eventually limit the internal mechanisms for deep convection, 
allowing for strong dry air entrainment that inhibits cyclone development. 

The Usagi system became detached from the baroclinic system and 
continued to move southwest between the two subtropical highs before finally 
moving west under the subtropical ridge and then northwest towards Japan. This 
system quickly strengthened as it approached the Mariana Islands, and the 
JTWC upgraded it to TS Usagi on 0000 UTC 29 July eventually becoming 
Typhoon Usagi 24 hours later. The JTWC advisories did not mention the 
possibility of SI from a decaying baroclinic system or the interaction with the 
remnants of Cosme. The first mention by JTWC of the seedling system was on 
27 July as a system was already undergoing SI. We believe that the prediction 
of this tropical cyclogenesis could have been anticipated farther in advance, if 
forecasters had more knowledge on Sis in the WNP. This case study, along 
with future work on others in the WNP will lead to a better understanding of this 
challenging area of tropical cyclone forecasting. 

The analysis suggest that the Usagi formation was the result of a SI from 
a strong baroclinic system in the mid-to-lower troposphere as it undergoes decay 
as it moved into an area of higher sea surface temperatures, weak low-level 
vertical wind shear, and persistent area of convection. As the PV anomalies 
were stretched and detached from the baroclinic source region, they were 
wrapped around a strong tropospheric anticyclone in the subtropics. This 
constitutes a different type of baroclinic initiation process than has been 
previously identified in Atlantic cyclone formation events associated with TT, 
which are induced by upper-level troughs. 
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The marsupial paradigm provides a convenient and self-consistent 
framework for the coupled dynamical/thermodynamical processes that are 
essential to tropical cyclogenesis. The study has highlighted the role curvature 
vorticity and moisture in the formation of Usagi through the analysis of the OW 
and SF parameters as done in DMW, 2008. We emphasized the formation of a 
closed gyre at the tip of the PV streamer from the baroclinic system, and 
confirmed that translated streamlines can be used as a good approximation of 
parcel trajectories to study the Lagrangian nature of the storm formation.. 

B. FURTHER STUDIES 

Further research is recommended on SI cases to better understand this 
type of tropical cyclone formation in the WNP, and thus contribute improved 
tropical cyclogenesis forecast. 

Future studies may include model studies that execute “vortex surgery” 
that remove Cosme and to examine whether or not Usagi would have eventually 
developed even in the absence of Cosme. This approach would allow a 
quantitative determination of the role of this circulation in the baroclinic system 
and the development of Usagi. 

It is hoped that this study will aid Tropical Cyclone Structure’ 08 (TCS’08) 
and THORPEX (THe Observing system Research and Predictability Experiment) 
Pacific-Asia Regional Campaign (PARC) in identifying the role of remnants of 
tropical cyclones in influencing the downstream flow. Do other storms 
completing SI like Usagi have similar anomaly patterns? Can these anomaly 
patterns be determined objectively and in a timely manner to be of use to the 
forecaster? These questions will be difficult to answer until a more robust data 
set is developed. 

Further investigation and more qualitative documentation are desired to 
increase JTWC forecast capability and such synoptic patterns or anomalies 
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favorable for tropical cyclone formation and to provide guidance on additional or 
mesoscale factors that may enhance or inhibit tropical cyclone development in 
the WNP. 
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APPENDIX A. OW PLOTS 
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Streamlines of horizontal flow at 600 hPa for the genesis sequence of the non¬ 
developer from 1800 UTC 23 July to 1200 UTC 24 July. The shadings indicate 
high values of OW parameter (units: 10' 10 s' 2 ) as defined in Eq. (1.1). This 
quantity, like vorticity, is invariant with respect to translation, therefore identical in 
translating and moving frames. Red areas are strong values and blue values. 
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Streamlines of horizontal flow at 600 hPa for the genesis sequence of the non¬ 
developer from 1800 UTC 23 July to 1200 UTC 24 July. The shadings indicate 
high values of OW parameter (units: 10' 10 s' 2 ) as defined in Eq. (1.1). 
This quantity, like vorticity, is invariant with respect to translation, therefore 
identical in translating and moving frames. Red areas are strong values and blue 
values areas are weak. 
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Streamlines of horizontal flow at 600 hPa for the genesis sequence of the non¬ 
developer from 1800 UTC 23 July to 1200 UTC 24 July. The shadings indicate 
high values of OW parameter (units: 10' 10 s' 2 ) as defined in Eq. (1.1). This 
quantity, like vorticity, is invariant with respect to translation, therefore identical in 
translating and moving frames. Red areas are strong values and blue values 
areas are weak of non-developer from 1800 UTC 22 July - 1200 UTC 25 July. 
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Streamlines of horizontal flow at 600 hPa for the genesis sequence of the non¬ 
developer from 1800 UTC 23 July to 1200 UTC 24 July. The shadings indicate 
high values of OW parameter (units: 10' 10 s' 2 ) as defined in Eq. (1.1). This 
quantity, like vorticity, is invariant with respect to translation, therefore identical in 
translating and moving frames. Red areas are strong values and blue values 
areas are weak. 
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APPENDIX B. RELATIVE VORTICITY PLOTS 
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Streamlines of horizontal (rotational+divergent) flow and Zeta at 600 hPa of the 
non-developer from 1800 UTC 22 July to 1200 UTC 25 July. Shading indicates 
relative vorticity (units: 10' 5 s' 1 ). 
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Streamlines of horizontal (rotational+divergent) flow and Zeta at 600 hPa of the 
non-developer from 1800 UTC 25 July to 1200 UTC 28 July. Shading indicates 
relative vorticity (units: 10' 5 s' 1 ). 
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Streamlines of horizontal (rotational+divergent) flow and Zeta at 600 hPa of the 
non-developer from 1800 UTC 22 July to 1200 UTC 25 July. Shading indicates 
relative vorticity (units: 10' 5 s _1 1). 
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Streamlines of horizontal (rotational+divergent) flow and Zeta at 850 hPa of the 
Usagi system from 1800 UTC 25 July to 1200 UTC 28 July. Shading indicates 
relative vorticity (units: 10' 5 s' 1 ). 
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APPENDIX C. SF PLOTS 
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Streamlines of horizontal flow 600 hPa for the genesis sequence of the non¬ 
developer, with high values of saturation fraction indicated by shading (unit: 
percent) from 1800 UTC 22 July to 1200 UTC 25 July. Red areas are strong 
values and white areas are weak values of SF. 
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Streamlines of horizontal flow 600 hPa for the genesis sequence of the Usagi 
system, with high values of saturation fraction indicated by shading (unit: 
percent) from 1800 UTC 25 July to 1200 UTC 28 July. Red areas are strong 
values and white areas are weak values of SF. 


92 













































































APPENDIX D. 


TRMM PLOTS 
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Streamlines of horizontal flow 850 hPa for the genesis sequence of the non¬ 
developer, with high values of TRMM 3-hour accumulated precipitation indicated 
by shading from 1800 UTC 22 July to 1200 UTC 25 July. Red areas are strong 
values and white areas are weak values of precipitation. 
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Streamlines of horizontal flow 850 hPa for the genesis sequence of the Usagi 
system, with high values of TRMM 3-hour accumulated precipitation indicated by 
shading from 1800 UTC 25 July to 1200 UTC 28 July. Red areas are strong 
values and white areas are weak values of precipitation. 
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